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. Abstract 

' We present a pedagogical description of the inversion of Gamow's tunnelling formula 

çSj ■ and we compare it with the corresponding classical problem. We also discuss the issue of 

$h . uniqueness in the solution and the result is compared with that obtained by the method 

Q" 1 ! of Gel'fand and Levitan. We hope that the article will be a valuable source to students 

■ who have studied classical mechanics and have some familiarity with quantum mechanics. 

o : 

'• 1 Introduction 

> : 

Eugen Merzbacher has commented that pQ "among all the successes of quantum mechanics 
çnj ■ as it evolved in the third decade of the 20th century, none was more impressive than the 

understanding of the tunnel effect — the penetration of matter waves and the transmission of 
ly-j ! partides through a high potential barrier." The tunnel effect provi ded a straightforward and 



remarkable explanation of the radioactive a-decay of nuclei. George Gamow was one of the- 
although not the sole- protagonists in the discovery of the theory of a-decay PJ Ej and the bàsic 
formula, equation (JT|). that underlies tunnelling through a potential barrier is often referred to 



as — perhaps unjustly — Gamow's penetrability factor (see, for example, p. 62 of |3j). 

Leaving aside Gamow's mischievous account of history, this article will discuss how knowl- 
qh| edge of the tunnelling behavior of a potential can be used to determine the potential itself. The 
quantum mechanical problem will also be contrasted with the classical version in order to gain 
further insight. 



2 Classical vs Quantum Problem 

2.1 The Classical Problem 

The classical, one-dimensional inverse scattering problem is probably best exemplified by the 
two systems depicted in flgure 

To the left of flgure ^ is a particle oscillating in an attractive potential with single minimum, 
set to zero and occurring at the origin. Conservation of energy gives the period, 



T(E) 




dx 



l( E) y/E-U{x) ' 



where X\(E) and X2(E) are the turning points for the energy E. The inverse problem is to 
determine the form of the potential, U(x), given the period as a function of energy, T(E). The 
solution is well known and may be found in many Standard texts on classical mechanics (see, 
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Figure 1: A particle in a potential well (attractive force field) performs oscillations; a particle 
on a potential barrier (repulsive force field) will either overcome the potential if it has enough 
energy (E > Uq) or will reflect back if it does not have enough energy E < Uq. 

for example [3]). Treating x as a function of U rather than U a function of x, the above may 
be turned into Abel's integral equation. Since U(x) is not one-to-one this requires splitting 
its domain at the origin and defining the two functions X\(U) and X2{U) as per figure ^ The 
result is 

(n\ fTA 1 í u T{E)àE 

x 2 {U) - xi(U) 



irVïmJo VU - E 

We find that the solution cannot be determined uniquely unless the additional, assumption 
that the potential is even, is introduced. In this case 

x{ü) 1 r 6 T(E)dE 



27T V / 2m Jo y/ü _ e ' 

where we have denoted the unique, even solution by U(x). 

To the right of figure ^ is a particle incident on a potential barrier that is confined to the 
interval [0, L\. The problem of determining the potential given the time of traversal of the 
potential as a function of energy has been solved by Lazenby and Griffiths [Sj. The forward 
and backward scattering data are defined respectively as, 

T{E) = Mi" , àx , ÜE>U , 
{ ' V 2 J y/E - U{x) 

R(E) = M r m ^ =tL i{E<U^ 
V 1 V 2 J y/E - U(x) ~ 

where X\(E) is the left turning point. The former applies when the particle has energy exceeding 
Uo, the maximum of the potential, and gives the time required for the particle to traverse the 
potential. The latter applies when E < Uq and gives the time required for the particle to reach 
the turning point x±(E) or half the time taken for the particle to return to the origin. 

The solutions to the barrier equations are quite similar to that of the previous system. The 
most important feature is that a class of potentials are obtained. There is, however, a unique 
solution with the property that it increases monotonically over the interval [0, L] (and drops 
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discontinuously to zero at L). Lazenby and Griffiths call this the canonical potential and use it 
to represent the class of solutions. For instance the inversion of the backward scattering data 
is given by 



where U is the canonical potential. 

In their paper, Lazenby and Griffiths remark that it is curious that the above solutions are 
not determined uniquely whereas in the quantum mechanical analogue the solution is unique. It 
is further remarked that "given the transmission coeficient T (the probability that the particle 
will surpass the barrier) as a function of energy E, (the potential) may be recovered by the 
method of Gel'fand and Levitan". As will be discussed, this statement is not accurate as it 
stands. The transmission coeficient alone is not suficient to determine the potential. The 
transmission amplitude however is a complex function and carries more information. All this 
will be clarified in the following. 

2.2 The Quantum Mechanical Problem 

The quantum mechanical problem is depicted in figure El 



Figure 2: Quantum mechanical scattering of a right-moving particle that approaches a potential 
from left. It is assumed that the potential vanishes at large distances and, therefore, bound 
states appear only if there is a region with U < and appear only for negative energies. 

A particle is incident from the left, say, on a potential which limits to a constant as x — > ±oo, 
which we will set to zero. If the potential approaches zero (a constant) rapidly enough, the 
asymptotic form of the wavefunction for x —>■ ±oo are plane waves: 





U 




—ikx 



X 



X 



— OO 



where the energy of the particle E > and 
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If there èxits an interval over which U(x) < 0, then there is a discrete spectrum E n correspond- 
ing to bound states: 



ip(x) 



71. C- J 



-oo 



d n e KnX , x —> +00 , 



where E„ < and 



-2mE n 



h 2 



The scattering data for the inverse problem is comprised of the asymptotic coeficients b(k) 
and c n as well as the discrete eigenvalues K n . The potential is then uniquely constructed by 
way of the method of Gel'fand and Levitan jü] outlined in figure EU below. 



b(k), (c„, K n ) 
ï 

I 

ÜT(x, z) + F(x + 2) + f+°° + z)dy = 



Figure 3: The Gel'fand and Levitan method in a nutshell. 

The scattering data is used to construct the auxiliary function F(X). The auxiliary function 
is then used to define a linear Fredholm integral equation, called the Marchenko equation, 
for a second auxiliary function K(x,z) (third entry). The potential is then determined by 
taking the directional derivative along the line z — x. 



3 Gamow's Formula and its Inversion 
3.1 Gamow's Formula 

The applicable situation is illustrated in figure |U A particle is incident on a potential barrier 
with a single maximum, Uq (at x — 0, say), with energy E less than Uq. It is known that in 
quantum mechanics there can be a finite probability for the particle to exceed the potential 
despite the fact that E < U . 

Gamow's tunnelling formula gives a good approximation to the transmission coefficient, 
T(E), giving the probability for the particle to surpass the barrier. 

/ 2 r X2{E) , \ 

T(E) = exp -- / y/2m(U(x) - E) dx . (1) 

V h Jx 1 {E) J 

This formula can be proved easily by considering the barrier as an infinite sum of infinitely thin 
rectangular barriers (p. 219 of [Zj). However, this method, although it provides the correct 
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Figure 4: A particle incident on a potential barrier has a finite probability to overcome the 
barrier even if its energy E is below the maximum Uq of the barrier. 

result, is mathematically inconsistent. A mathematically sound proof can be given using the 
JWKB approximation (p. 507 of jZj). 

In what follows, T(E) will play the role analogous to the classical scattering data. 

3.2 Inversion of Gamow's Formula 

Let us take up the task of inverting Gamow's formula. By differentiating and, once again, 
rewriting in terms of the inverse functions x±(U) and X2(U) (we have again split the domain of 
U(x) at the origin), we find 



h i dT r X2(E) dx r Uo /dxi _ dx 2 \ du 

~^T\E)dE 1 J xm ~ J E \dÜ ~ dÜJ y/JT^E 

The above equation is nearly in the form of Abel's equation 



(2) 



Jo y/E^ÏJ 



However, it differs in that the position of the parameter and the variable have been switched in 
the root, and in that the limits of the integral are from the variable to a constant rather than 
zero to the variable: 



' e y/Ü^E 

Consequently, the preferred approach of applying the Laplace transform to the equation and 
making use of the convolution theorem fails. However, Abel's equation can be solved by com- 
position with a kernel jl]. We can, in fact, apply this process with some modification. We 
divide both sides by y/E — a, where < a < Uo, and integrate with respect to E from a to 

h f Uo dT/dE lrn f Uo dE f Uo fdxx dx 2 \ dU 

-dE 



2mJ a T(E)y/E-a J a y/E~^~a J e \dU dU ) y/\J - E 



'd Xl dx 2 \ du f U 



dU dUJ J a y/{U -E){E-a) 
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E = U 




Figure 5: Our integral is defined on the triangular domain T>. 



where in the second line we have changed the order of integration (see figure 
In the appendix we show that 



Therefore 



and finally: 



v 



dE 



y/(U-E)(E-a) 



h 



'2m 



U 



dT/dE 



T(E)VW 



-AE 



7T 



Uo 



dx\ dx- 



dU dU 



dU 



x 1 (U)-x 2 (U) 



h 



7rv2m 



u 



dT(E)/dE 



dE 



(3) 



T(E)VE- U 

where we have used the fact that xi(Uo) — 22(^0) = 0. We see that the solution is quite similar 
in form to the classical result, and again the solution is not unique. Rather, we have obtained 
a family of potentials which all result in the same transmission coefncient. 



3.3 An Example: Cold Emission 

Although the photoelectric effect has become the most popular example of emission of electrons 
by metals, easily recognized and studied by students, it is not the only phenomenon in which 
electrons are emitted by metals. Electrons can be emitted by metals at room temperature by 
the application of an external elèctric field 8. To contrast with the emission of electrons when 
a metal is heated, this phenomenon is termed cold emission. It was first explained by Fowler 
and Nordheim 

When an external field is applied, an electron in the metal sees a potential 

U(x) = U - e£ x , (4) 

where x is the distance from the wall of the metal. This description ignores the fact that a 
positive image charge will appear at the surface of the metal as the electron is removed, and 
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thus an additional Coulomb attraction will be established. For an electron of energy E we find, 
using Gamow's formula (JTJ), 



T(E) 



with 



-a([/ -E) 3/2 



4V2m 



(5) 



3e£h 



Equation (jSJ) is known as the Fowler-Nordheim equation. The quantity Uq — E is known as 
the work function. 



7\ U(x) = Uq — e£ x 




Figure 6: Model potential for cold emission. E is typically taken as the fermi energy. W is the 
work function of the metal. 



We shall now assume that the Fowler-Nordheim equation is known — say from experimental 
data. Can we find the potential that reproduces it? 
According to our formula (J3J) 

This is an elementary integral and we compute it in the appendix. The result is 

x 1 (U)-x 2 (U) = ~{U -U) . 
et 

The reader may believe that we have recovered the potential (j3J). Unfortunately, this is not the 
case as any two functions X\(U) and X2(U) that differ by the above amount are solutions of the 
inverse problem. The the cold emission potential is recovered if we assume that Xi(U) = 0. 



3.4 The Issue of Uniqueness 

The reader, at this point, may believe that the apparent conflict of our result with that which 
would have obtained by the method of Gel'fand-Levitan is due to the approximations used to 
produce formula However, a moment's reflection will reveal that this cannot be so. The 
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mathematical statement of the problem is independent of the underlying physics which can be 
left aside. 

The answer to this puzzle is quite simple. The method of Gel'fand and Levitan makes use 
of the amplitude b(k) (in the situation were considering there is no bound spectrum) which is a 
complex quantity. However, in our case we make use of T which is a real quantity — T = |a(/í;)| 2 
with |a(A;)| 2 + |fr(&)| 2 = 1. Thus, we have lost information about regarding phase. 

Given the transmission coefficient T(E), b(k) can be anything of the form 



where f(k) is a real-valued function. Each distinct potential among the family of our solutions, 
corresponds to a different choice of f(k). 

4 Discussion and Conclusion 

We have found that, in quantum mechanics, T(E), the probability for transmission and the 
analog of the classical scattering data, does not uniquely determine the potential, just as it is 
in classical mechanics. However, quantum mechanics does afford us an additional set of data, 
the phase difference f(k), which corresponds to measurements of time delay (p.138 of jS]). It 
is only with both these sets of data that we can uniquely determine the potential. 

The approximate nature of Gamow's formula is an irrelevant feature for the problem we 
have studied. However, a different kind of question may be asked which makes this feature 
relevant: Assuming that the potential barrier is even and a unique solution U(x) can be found, 
what is the error in determining the potential? That is, how close is the solution U(x) to the 
real potential, U rea i(x), that gave the experimental data T(E)1 This question remains open. 

Finally, toward proving our result, we have succeeded in solving a modified version of Abel's 
equation: Given the integral equation, 



where f(E) is a known function and 4>(U) is an unknown function, we have shown that the 
solution is given by 
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b(k) = y/l-T(E(k)) e im , 
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Appendix 

The integrals 



dE 



^((3-E)(E-a) 



, J 



dE 



which have appeared in our article are elementary, but their calculation is somewhat lengthy. 
We present their calculation here. 
Introducing the substitution 



tr 



13 -E 
E-a 



in J, we can rewrite it as 

J = -2((3-a) 
Howcver, 



ir 



(1 +u 2 ) 2 



du = -2(p- a) 



1 + u 2 



du 



(1 +u 2 ) 2 



du 



1 , u 
— tan — 
A A 



A 2 + u 2 



Therefore, 



and, finally, 



J 



2(j3-à) 



Q tan" 1 






fi 








A=l 



J = ^{E-a){(3-E) -{(3-a) tan^ 1 J 

V L·i — oi 



It is immediate that 



E — a 2 



Using Leibnitz's rule for the differentiation of integrals depending on a parameter, we can also 
easily obtain 

dE dl 

2— = 7T . 



- WÉ - a) ~d(3 
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